Surface effects in color superconducting strange-quark matter 
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Surface efTects in strange-quark matter play an important role for certain observables which have 
been proposed in order to identify strange stars, and color superconductivity can strongly modify 
these effects. We study the surface of color superconducting strange-quark matter by solving the 
Hartree-Fock-Bogoliubov equations for finite systems ( "strangelets" ) within the MIT bag model, 
supplemented with a pairing interaction. Due to the bag-model boundary condition, the strange- 
quark density is suppressed at the surface. This leads to a positive surface charge, concentrated in a 
layer of ~ 1 fm below the surface, even in the color-flavor locked (CFL) phase. However, since in the 
CFL phase all quarks are paired, this positive charge is compensated by a negative charge, which 
turns out to be situated in a layer of a few tens of fm below the surface, and the total charge of 
CFL strangelets is zero. We also study the surface and curvature contributions to the total energy. 
Due to the strong pairing, the energy as a function of the mass number is very well reproduced by 
a liquid-drop type formula with curvature term. 

PACS numbers: 21.65.Qr,12.39.Ba,26.60.-c 
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I. INTRODUCTION 

From rather general arguments it is expected that at 
low temperatures and high densities quark matter is in 
a color superconducting state More recently 0, Q 
it has been suggested that the diquark pairing gaps for 
quark matter at densities of several times nuclear mat- 
ter saturation density could be of the order of ~ 100 
MeV. Since this could have important phcnomenologi- 
cal consequences in particular for the interior of compact 
stars, this has triggered much work on color supercon- 
ductivity in dense quark matter (for reviews, see, e.g., 
Ref. [1]). These investigations of the QCD phase diagram 
have revealed a very rich phase structure with many dif- 
ferent possible pairing patterns, depending on external 
conditions such as, for instance, electrical neutrality or 
quark masses. The largest diquark pairing gaps arise 
from scalar condensates, leading either to the two-flavor 
color superconducting_(2SC) phase or to the color-flavor- 
locked (CFL) phase [a, Q. The latter pairing pattern 
involves strange (s) quarks, in addition to the two light 
quark flavors, up (u) and down (d). 

If color superconducting quark matter exists in nature, 
the most likely place to find it is the interior of com- 
pact stars because matter is compressed there to densi- 
ties much higher than nuclear matter saturation density. 
However, it has been argued that strange-quark matter 
(SQM) might be absolutely stable [7,]. Under this hy- 
pothesis, even pure strange stars should exist P, i.e., 
stars entirely composed of SQM. Also small lumps of 
SQM, called "strangelets," might be stable. Because of 
their low charge to baryon number ratio Z/A, strangelets 
have been proposed to populate ultra- high energy cosmic 
rays 

In SQM without pairing, the density of strange quarks 
is supposed to be smaller than that of light quarks be- 



cause of their higher mass. Consequently, SQM and 
strangelets are positively charged and the charge neutral- 
ity of strange stars has to be achieved via the presence 
of electrons. At the surface an atmosphere of electrons 
forms P which can potentially be detected [13, [HI via 
the emission of electron-positron pairs from an extremely 
strong electric field at the surface. 

Recently another possible picture of the surface of a 
strange star has been proposed there could be a 

"crust" composed of strangelets immersed in an electron 
gas. Similar to an ordinary neutron star, there could be 
an interface between the crust and the interior in form 
of the famous "pasta phases." Within this scenario the 
electric field at the surface would be strongly reduced. 
Obviously, surface effects for the strangelets play an im- 
portant role for the description of this scenario. For in- 
stance, there is a critical surface tension deciding whether 
a homogeneous phase or the droplet phase is favored [l^ . 
Another question for which surface effects should be con- 
sidered is the formation of a strange star in a supernova 
explosion. Before the explosion the original star contains 
hadronic matter. During the formation of the star, nucle- 
ation of strangelets sets in, leading then to a conversion 
of the entire star to SQM. For the nucleation process the 
properties of small strangelets are important. 

Pairing tends to reduce the differences in density of dif- 
ferent quark species. For bulk quark matter in the CFL 
state, requiring color neutrality, all quarks are paired. 
The densities are thus equal and CFL quark matter is 
electrically neutral on its own, i.e., without any elec- 
trons [14]. This would suggest dramatic changes in the 
properties of strangelets and SQM inside compact stars. 
For instance, the electrosphere at the surface of a strange 
star could completely disappear. But, the presence of 
the surface can modify this picture since it can lead to 
a non-zero surface charge which remains even for large 
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objects. For example, the boundary condition of the 
MIT bag model suppresses the density of the massive 
strange quarks at the surface, resulting in a positive sur- 
face charge [l5|. Within this scenario, the total charge 
of a strangelet, following roughly Z k, O.SA^/"^, is drasti- 
cally reduced with respect to "normal" strangelets. For 
strange stars, this requires the presence of electrons 

PI- 

However, pairing has not been treated self-consistently in 
previous work (see, e.g., Ref. [IBl)- In this paper we will 
therefore reinvestigate finite-size strangelets with pair- 
ing by considering quark matter in a color superconduct- 
ing spherical bag, solving the Hartree-Fock-Bogoliubov 
(HFB) equations. We will show, in particular, that there 
exist CFL type solutions where all quarks are paired and 
the total charge of the strangelet strictly vanishes. 

The outline of the paper is as follows. In Section [Til we 
will present our model for treating color superconduct- 
ing quark matter in a finite volume. In Section IIIII we 
will show numerical results. In Section IIII Al we discuss 
the possibility of qualitatively different configurations. In 
Section IIII Bl we concentrate on the charge-density dis- 
tributions of the CFL like solutions. In Section IlII CI we 
discuss a liquid-drop like mass formua for the CFL-like 
solutions and calculate the surface tension. Finally, in 
Section ITVl we will summarize our results. 



II. MODEL 

A. Lagrangian 

Since it is not possible to describe strangelets or SQM 
with a surface from first principles (QCD), we will use 
a quark model which allows to describe finite-size ob- 
jects. For this purpose we will use here the MIT bag 
model The idea of this model is that confinement 

can be simulated by the existence of a "bag" which con- 
sists of a "hole" in the non-perturbative QCD vacuum. 
Inside this "bag" , the vacuum is supposed to be pertur- 
bative, i.e., inside the bag the interactions of the quarks 
can be treated perturbatively. To create this "hole" in 
the non-perturbative QCD vacuum, an energy per vol- 
ume, i?, is necessary. In the present work we will consider 
a static spherical bag with radius R. On the surface of 
the bag, the quark field '0 h^^s to satisfy an appropriate 
boundary condition. In the simplest version of the MIT 
bag model, the boundary condition reads 

- ie^ • 7-0 = -0^=^ • (1) 

which ensures that there is no particle flux across the 
surface. By r = |r| we denote the radial coordinate, 
measured from the center of the bag, and = r/r is the 
radial unit vector. The boundary condition ([T]) leads to a 
suppression of the wave functions of massive particles at 
the surface. This means that the strange-quark density 
will a priori be suppressed at the surface with respect to 
the light quark densities. 



The MIT bag model can be expressed in terms of a 
Lagrangian density as follows 

Cbag = miY'd^ ~ m)^ - B]0{R - r) 

- iv^V-^li? - r) , (2) 

where m is the matrix of quark masses. Due to the sec- 
ond term, the boundary condition ([1]) follows immedi- 
ately from the Euler-Lagrange equation for the quark 
field [3. 

In order to include pairing, we will supplement the bag 
model with a pairing interaction. In principle, perturba- 
tive one-giuon exchange generates an attractive paring 
interaction in certain channels, in particular in the scalar 
color antitriplet channel. For simplicity, we will use here 
a four-point pairing interaction acting only in this dom- 
inant channel. The corresponding Lagrangian reads (see 
any of the standard review articles on color superconduc- 
tivity [4]) 

>Cp„r = ^ X! (V'«75TAAA'C'V5^)(0^Ci75TAAA'V') , (3) 
A,A' 

where _ff is a dimensionful coupling constant, C de- 
notes the charge conjugation matrix, and ta,Xa' rep- 
resent S'?7(3) matrices in flavor and color space, respec- 
tively. We follow the convention that capital letters A, A' 
indicate that we are restricting ta and Xa' to be an- 
tisymmetric, i.e., in terms of the Gell-Mann matrices, 
A^'G {2, 5,7}. 

In addition to the strong interaction, the quarks will 
exhibit electromagnetic interactions which, due to their 
long range, become particularly important for large ob- 
jects. The corresponding Lagrangian reads 

Ce.m. = '-\f^.F^''' - eijQA^Yi', (4) 

where F^^^ = d^j^A^—d^Af^ and denote respectively the 
electromagnetic field strength tensor and four-potential, 
and Q is the matrix of quark charges in units of e, Q„ = 
2/3, Qd^Qs^ -1/3. 

It would be in the spirit of the bag model to include 
also the giuon exchange in a perturbative way, i.e., in the 
same way as the photon. However, this goes beyond the 
scope of the present paper and will be postponed to a 
future study. 

B. Solution in the framework of HFB theory 

The model described by £ = Chag + C,pair + ^e.m. will 
be treated in the framework of HFB theory. By mini- 
mizing the energy in mean field approximation (for more 
details see Appendix |B] and Ref. [19|] where the "Dirac- 
Hartree-Bogoliubov" approximation was developped for 
finite nuclei), one obtains the following HFB equations: 

(h ^\( C/„(r) \ _f f/.(r) \ . . 
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The single-particle Hamiltonian 

h = -ia ■ V + ?7i7° + S - /i (6) 

includes besides the free Dirac Hamiltonian the quark 
self-energy E (in our case due to Coulomb interaction) 
and the matrix of chemical potentials fi which depend 
on flavor / € {u,d,s} and color c G {r,g,b} (we will 
denote the three colors by red, green, and blue). A de- 
notes the pairing field (gap). The spinors Ua and Va 
describe the particle- and hole-like components of the 
quark fields, respectively [see Eq. (|B3|) ]. where a is a 
multi-index containing all quantum numbers character- 
izing a single-particle state (see Appendix [X)) . In writ- 
ing Eq. ([5]), we implicitly assumed that the pairing field 
A can be chosen real, which is the case for the pairing 
pattern we consider, and that the self-energy E is local, 
which is the case since we neglect the exchange (Fock) 
term (see below). 

The pairing field A and the self-energy E depend them- 
selves on the wave functions U and V, such that we have 
to solve a self-consistency problem. To be specific, the 
pairing field A depends on the diquark condensates 

SAA'{x) = - {TfT{x)TAXA''ip{x)) , (7) 

where tpx denotes the time-reversed conjugate of V', 

V'T = IsC^'^ ■ (8) 

The diquark condensates can be expressed in terms of 
the U and V functions as 

SAA'{r) = - J2 Vf3ir)TAXA'Upir) (9) 

(since we are dealing with a static problem, the con- 
densates do not depend on time, and due to spherical 
symmetry, they depend only on the radial coordinate r) . 
We will limit our investigations here to diagonal conden- 
sates, i.e., only condensates with A = A' are non-zero"'^. 
In uniform infinite matter and for an exact SU (3) fiavor 
symmetry, the CFL phase is characterized by nonzero 
values S22 — S55 — sjt, whereas the 2SC state has only 
S22 7^ 0. The relation between the condensates saa and 
the pairing field A reads 

A(r)= ^A{r)TAXA, (10) 

A=2,5J 

AA{r) = 2HsAA{r). (11) 

In practice, the expression ^ is divergent and it is nec- 
essary to introduce a cutoff in order to obtain a finite 



^ In uniform infinite matter it can be shown that for the ener- 
getically favored solution the arbitrary orientation in color can 
be chosen in such a way that only the diagonal condensates with 
A = A' are non-zero. 



result. Since in a finite system the levels are discrete, a 
sharp cutoff would generate discontinuities as a function 
of the system's size. We therefore introduce a smooth 
cutoff function f{p/ A) (see Appendix [Cl for details). An- 
other practical problem arises from antiparticle contri- 
butions. However, since the chemical potentials /i/c are 
large and positive and pairing involves mostly the states 
near the Fermi surface, we assume that the antiparticle 
contributions are not important and can be neglected. 
We checked this approximation (analogous to the "no-sea 
approximation" in nuclear physics 19]) in infinite mat- 
ter and found that the effect of antiparticle states can be 
absorbed in a readjustment of the coupling constant by 
- 20 %. 

For the normal self-energy E we employ the Hartree 
approximation, i.e., we neglect the Coulomb exchange 
(Fock) term as well as exchange contributions from the 
magnetic field. We also disregard the contribution of 
^pair to the normal self-energy. Hence, the self-energy is 
simply proportional to the static Coulomb potential 

E(r) = eQ^(r)7° . (12) 

The Coulomb potential is related to the quark densities 
by 

Mr)-ejd^r'^^, (13) 

where 

Pc/.(r)=^Q/P/(r) (14) 
/ 

is the charge density (divided by e), being the number 
density of quarks of fiavor /. As it was the case for the 
diquark condensates, the quark number densities can be 
expressed in terms of the U and V functions. Denoting 
by (3 all single-particle quantum numbers except flavor, 
we can write the number density of quarks of fiavor / as 

Let us now summarize the procedure how the HFB 
equations are solved. We start with an initial guess for 
the pairing fields AA{r) and for the Coulomb potential 
Ao{r). Then we solve the eigenvalue problem (O in order 
to find the U and V functions. From these functions 
the diquark condensates SAAif) and the quark densities 
Pf{r) are computed according to Eqs. ([9]) and (fT5|) . which 
are then used to update the pairing fields A^ and the 
Coulomb field Aq according to Eqs. pT|) and p3|l . These 
steps are iterated until convergence (i.e., self-consistency) 
is reached. 

The crucial difference to the BCS formalism in ho- 
mogeneous infinite matter is that in our case the wave 
functions adapt themselves to the pairing field and to 
the Coulomb potential, whereas in the case of homo- 
geneous infinite matter the wave functions always stay 
plane waves, and the U and V factors are just coefficients 
multiplying them. 
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Determination of chemical potentials and bag 
radius 



In Section IIIBI we described how the HFB equations 
are solved for given values of the chemical potentials /i /c 
and of the bag radius R. However, in reality, only one 
quantity is given, namely the baryon number A. Even 
the fractions of different quark flavors cannot be fixed, 
unless one allows for /3 unstable strangelets. Let us now 
describe how we determine the chemical potentials fife 
and the bag radius R for given baryon number A. 

The first step consists in fixing the quark numbers, 
Nfc, for each flavor / and color c, and to adjust the che- 
mical potentials fi fc in order to obtain these quark num- 
bers. Before we address the question how the nine quark 
numbers Nfc are determined, let us discuss the issue of 
the bag radius R. Until now, the radius was imposed 
from outside, but in reality the system will choose its 
radius such that it minimizes its total energy for given 
quark numbers Nfc- Within the bag model, the total 
energy is given by 



E = Ea + BV . 



(16) 



where V — 47ri?^/3 is the volume of the bag. By Eq we 
denote the energy of the quarks inside the bag, includ- 
ing the interaction energy, which in our case comes from 
pairing and Coulomb interactions. It can be obtained 
from the solution of the HFB equations as follows [T9| : 

/ '^''^ E {ul{v){ep + fi)Uf3{v) 



+ 2 



[/Mr)A(r)V}3(r)-[/T(r)eQAo(r)(7^(r) . (17) 



Minimizing the total energy E is of course completely 
equivalent to saying that the quark pressure in the bag 
is counterbalanced by the bag pressure B, i.e., 



dEq 
dV 



(18) 



N 



This equation determines the radius of the strangelet for 
given bag pressure B, interaction strength H and quark 
numbers Nfc- In practice, however, we flnd it more con- 
venient to minimize E rather than solve Eq. ^TE\i . 

Let us now turn to the determination of the quark 
numbers. The nine quark numbers Nf^ cannot be chosen 
arbitrarily, but they have to fulfil certain requirements. 
Imposing the total baryon number A and color neutrality, 
i.e., equal numbers of quarks for each color, we have to 
satisfy the constraint 



A 



for all c . 



(19) 



unique values for the Nf^ it is necessary to impose /3 
stability, as we will describe now. 

In an infinite homogeneous system the condition for 
f3 equilibrium gives just a relation between the chemical 
potentials^ 



f^dc /^sc ^uc ~t- /ie 



for all c . 



(20) 



In a small system this is slightly different. First, even if 
there are electrons (i.e., if the strangelet is charged), they 
are not localized inside the strangelet, but they form a 
large cloud like in ordinary atoms and hence their che- 
mical potential fi^ is approximately equal to the electron 
mass and can be neglected. Without pairing, it has been 
estimated in Ref. [20| that this may be still the case for 
strangelets with charge Z < 1000, corresponding roughly 
to A < 10^. The second difference to bulk matter comes 
from the fact that, due to the discrete levels, particle 
numbers are discontinuous functions of the chemical po- 
tentials. The term (3 equilibrium should now be replaced 
by P stability, which means that the system does not 
gain energy by performing a (3 decay, inverse /? decay, or 
electron capture, i.e., transforming an up into a down or 
strange quark, or vice versa, accompanied by the corre- 
sponding leptons. 

To achieve /3 stability, we therefore compare the en- 
ergies of adjacent strangelets with the same total quark 
number per color, differing only in the number of up, 
down, and strange quarks, respectively, in order to 
find the configuration with the lowest total energy E. 
Of course, in the case of large particle numbers, the 
minimum-energy configuration fulfils approximately the 
condition ^U\j. 



D. Choice of the model parameters 

Besides the quark masses, which we take as to„ = 
rrid = and = 120 MeV, our model contains three 
parameters: the bag constant B, the coupling constant 
of the pairing interaction, H, and the cutoff A which is 
necessary to avoid the divergence of the gap equation 
see below Eq. (fTTj) . In fact, a change of the cutoff 
in reasonable limits can to very good approximation be 
compensated by a change of the coupling constant. We 
therefore choose rather arbitrarily A = 600 MeV and give 
instead of the dimensionful coupling constant H the di- 
mensionless combination HA^. So we are left with two 
parameters, B and HA" 

We can get an idea of the value of the bag pressure 
by looking at the stability of bulk quark matter. Non- 
strange quark matter should be energetically less favored 
than normal hadronic matter, whereas SQM should be 
stable if for some baryon number A > Ac strangelets 



Of course, these three equations are not sufficient for de- 
termining all the nine quark numbers. In order to get 



^ Here we assume that neutrinos are not trapped, i.e., they can 
freely leave the system 
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TABLE I: Values of the bag constants for different values of 
the coupling constant H, resulting in color and electrically 
neutral SQM with electrons in /3 equilibrium with an energy 
per baryon of E/A = 900 MeV. The corresponding baryon 
densities ps, electron densities pe, and pairing gaps in infinite 
matter are also displayed. 





(MeV) 


PB 

(fm-=^) 


(fm-^) 


Aa 
(MeV) 


As = A 
(MeV) 





152.03 


0.329 


7.3x10-'' 








1.5 


152.44 


0.339 


9.7x10-^5 


27.7 





1.75 


153.97 


0.367 





35.1 


34.5 


2 


156.26 


0.395 





50.6 


49.7 


2.25 


159.46 


0.427 





67.2 


66.0 


2.5 


163.46 


0.463 





84.6 


83.1 



0.5 




1 2 3 4 5 
r (fm) 

FIG. 1: Quark number density profiles of the strangelet A = 
108, Z = 24 in the case of vanishing pairing interaction (free 
quarks in a bag) and B^'^ = 152.03 MeV. 



become stable and consequently strange stars can exist. 
This means that we want the energy per baryon of SQM 
to be less than 931 MeV, the energy per baryon of the 
most stable nucleus, ^^Fe. On the other hand, the energy 
per baryon of non-strange quark matter should be larger 
than the nuclcon mass. Without interaction the win- 
dow for the values of the bag constant is then 148 MeV 
< JS^/^ < 157 MeV. These values change as a function 
of the interaction strength H. To better compare the 
results, we will readjust for each coupling strength the 
bag constant in order to get E/A — 900 MeV. The corre- 
sponding values are listed in Table [H together with other 
properties of infinite matter. Non-strange quark matter 
is unstable with these parameter values. Note that for 
the weakest non- vanishing coupling constant given in Ta- 
blelH SQM is in the 2SC phase and not in the CFL phase. 
For the larger coupling constants, the CFL phase is pre- 
ferred. Note that, due to the mass difference of light and 
strange quarks, the flavor SU{?>) symmetry is not exact 
and the gap A2 is different from A5 and A7. However, 
since the CFL phase is electrically neutral, and we have 
rriu = fTid — 0, the isospin SU{2) symmetry in the up- 
and down-quark sector is exact and therefore A5 — A7. 



III. RESULTS 

A. Different types of solutions 

We will first discuss the qualitatively different con- 
figurations we find. Let us start by discussing a small 
strangelet {A — 108, Z = 24) without any pairing 
interaction {HA^ = 0). The mass number has been 
chosen such that the minimum-energy configuration is 
a closed-shell configuration. The quark numbers and 
other relevant information are listed in Table |TT1 Due 
to the finite size of the bag, the energy per baryon 
(E/A = 932.5 MeV, including 1.0 MeV due to Coulomb) 
is much higher than that of color neutral infinite matter 



with flue = fJ.dc = fJi-sc^ {E/A = 899.5 MeV). This effect 
will be discussed in more detail in Section IIII CI The 
density profiles of light and strange quarks are shown in 
Fig. [1] As expected, due to the boundary condition, the 
strange-quark density is strongly suppressed at the sur- 
face, contrary to the densities of the light quarks. For 
comparison we mention that for the same value of the 
bag constant, the densities in color neutral infinite mat- 
ter with flue = Mdc = Msc are: Pu = Pd ^ 0.355 fm"^, 
Ps = 0.274 fm"^. We see that not only the strange-quark 
density, but also the densities of the light quarks are quite 
different from these values and depend strongly on r be- 
cause of the existence of discrete levels in the bag. Let us 
mention that, due to the Coulomb potential, the density 
profiles of up and down quarks are slightly different, but 
the difference is too small to be visible in Fig. [1] 

Now we switch on the pairing interaction. In the case 
of HA'^ = 1.5, SQM is in the 2SC phase, i.e., only up 
and down quarks of two colors (red and green in our no- 
tation) are paired. This is also true in a finite strangelet. 
Therefore it is clear that the strange-quark density pro- 
file remains the same as without pairing. The oscilla- 
tions of the densities of the light quarks, however, are 
much weaker now than in the case without pairing, since 
pairing washes out the occupation numbers. This can 
be seen in Fig. [51 In this 2SC-like solution, only one of 
the gaps, A2, is non-zero. Since A2 involves only the 
wave functions of up and down quarks, which are not 
suppressed at the surface, it extends up to the surface 
of the bag, as shown in Fig. [31 As a function of r, it 
is almost constant and quite close to the corresponding 
value in infinite matter with puc — fJ-dc — fJ-sc, which is 
A2 = 29.2 MeV. 

If we increase the coupling constant to HA^ — 1.75, we 



^ As discussed below Eq. I I20I I. it is more appropriate to compare a 
small strangelet with this kind of matter rather than electrically 
neutral matter with electrons in /3 equilibrium. 
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TABLE II; Parameters and properties of the strangelets discussed in Section Till Al B = bag constant, H — coupling constant 
of the pairing interaction, A = baryon number, Z = charge, Nfc = number of quarks of flavor / and color c, E/A = energy 
per baryon, R = radius of the bag, Aa(0) = value of the gap at r = 0. 



Bi/4 




A 


Z 


/ ]V„, ]V„g ]V„(, \ 

N^r Ndg 


E/A 


R 


A2(0) 


A5(0) 


A7(0) 


(MeV) 


(MeV) 




\ ]V„ iV„ JV,(, J 


(MeV) 


(fm) 


(MeV) 


(MeV) 


(MeV) 



152.03 108 24 (tttttt) 932.5 4.36 

V 20 20 20 / 



152.44 1.5 108 24 (tiiHi) 930.7 4.31 32.9 

V 20 20 20 / 

153.97 1.75 108 24 ft! III!) 934.0 4.24 49.5 

V 20 20 20 / 

153.97 1.75 108 10 f 39 1 934.8 4.21 41.6 24.9 24.9 

V 31 31 26 / 

153.97 1.75 108 f i i i ) 934.8 4.17 33.9 37.0 36.9 

V 35 35 38 / 
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FIG. 2: Quark number density profiles of the strangelet A — 
108, Z = 24 in the case of HA'^ = 1.5 and B^/" = 152.44 
MeV. 
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FIG. 4: Density profiles of the strangelet A = 108, Z = 10 in 
the case of HA'^ = 1.75 and B^^* = 153.97 MeV. 
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FIG. 3: Gap A2(r) of the strangelet A = 108, Z = 24 in the 
case of HA'^ = 1.5 and B^/* = 152.44 MeV. 



obtain three qualitatively different solutions which have 
comparable energies. The most stable one is still of the 
2SC type, although in infinite matter the CFL phase is 
preferred. In this case, the strangelet still has Z — 2A and 
the density profiles are almost identical to those shown 
in Fig. O The main difference is that now the value of 



the gap is larger. 

In the two other solutions, also strange quarks partici- 
pate in pairing (A5 w A7 7^ - note that A5 and A7 are 
not exactly equal because the isospin symmetry is broken 
by the Coulomb interaction). These two solutions have 
charge Z = 10 and Z = 0, respectively. Let us first dis- 
cuss the case Z = 10. In this case, there are a couple of 
up and down quarks which remain unpaired. The wave 
function of the unpaired level is mainly localized near the 
surface of the bag, as can be seen in Fig. [U where the 
density profiles are shown. In the inner part, the densi- 
ties of up, down, and strange quarks are almost equal, 
while near the surface, where the strange-quark density 
is suppressed due to the boundary condition, there is an 
excess of up and down quarks. This excess is due to 
the unpaired quarks. The fact that one level of up and 
down quarks (in the present case the 139/2 level, i.e., the 
lowest level with j = 9/2, k = — 5 in the notation of Ap- 
pendix]^ does not participate in pairing means that the 
occupation number of this level is equal to 1. At the same 
time, the corresponding level of the strange quarks has 
an occupation number equal to 0. In a certain sense this 
situation is analogous to the "breached pairing" phase 
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FIG. 5: Gaps Aa as functions of r for the strangelet A — 108, 
^ = 10 in the case of HA'^ = 1.75 and B^/* = 153.97 MeV. 
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FIG. 6: Density profiles of the strangelet A = 108, Z — in 
the case of HA'^ = 1.75 and B^^^ = 153.97 MeV. 

of infinite matter [2l[. The charge Z is equal to the de- 
generacy 2j + 1 of the unpaired level. The gaps as 
functions of r corresponding to this solution are displayed 
in Fig. O 

In the third solution, all quarks are paired. As a con- 
sequence, the numbers of up, down, and strange quarks 
are equal, and the total charge is Z — 0. This is analo- 
gous to the CFL phase in the infinite system. Since the 
strange-quark density is suppressed near the surface, but 
the number of strange quarks is equal to that of up and 
down quarks, it is clear that the strange-quark density 
must be larger than the up- and down-quark densities in 
some other part of the system. This is indeed the case, 
as can be seen in Fig. [SI We also see that the excess of 
the light-quark densities over the strange-quark density 
is reduced as compared with the case Z = 10 discussed 
above (cf. Fig. 0]). We will discuss the charge-density 
distribution in detail in Section fill Bl The gaps, shown 
in Fig. [71 are much closer to the gaps in infinite matter 
(cf. Table [J) than in the case Z = 10. 

For the larger values of the coupling constant we con- 
sidered (FA^ = 2,2.25,2.5), it is always the CFL-type 
solution {Z — 0) which has the lowest energy. We do not 
show any figures because in all these cases the results are 
analogous to those shown in Figs. [n|and[7| (just the va- 



lues of the gaps change, they are close to those given in 
Table [J for infinite matter). 

It should be mentioned that the fully paired solutions 
with Z — are very robust as soon as the coupling con- 
stant is sufficiently large, i.e., we find this type of so- 
lution for arbitrary numbers of quarks'*. This solution 
is in contrast to previous findings (see, e.g., Ref. JLS]), 
where it was supposed that the CFL matter should be 
neutral in the bulk with just a thin positively charged 
surface layer with an excess of up and down quarks be- 
cause of the boundary condition. In fact, this idea cor- 
responds roughly to our solution with unpaired up and 
down quarks near the surface. This solution is, however, 
very fragile and exists only for certain values of parame- 
ters and mass numbers, since it requires the existence of 
a suitable level of light and strange quarks near the re- 
spective Fermi surfaces which can serve as unpaired level. 



B. Charge density distribution 

We have seen in Section IIIl Al that in all cases ex- 
cept the 2SC phase, pairing drastically reduces the total 
charge Z. Because of surface effects, the local charge den- 
sity does, however, not vanish, even within the CFL-type 
solution which has Z — 0. Due to the suppression of the 
strange-quark wave function at the surface, a positively 
charged surface layer remains with an extension of ~ 1 
fm, as has already been pointed out in Ref. [l5l |. 

Within the configuration with some unpaired light 
quarks at the surface, the total charge of the strangelet 
results from this positive surface charge, the interior of 
the strangelet has almost zero charge density. The to- 
tal charge is here reduced compared with a strangelet 
without pairing, for example the A = 108 strangelet has 



* If the number of quarks is odd, it it impossible to pair all quarks 
and one or several state (s) should be "blocked" by the unpaired 
quark(s). At present, we have not included this effect in our 
calculation, and we restrict ourselves to even quark numbers 
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FIG. 8: Charge density profiles of tfie fully paired {HA'^ = 2) 
strangelets A =108, 3000, 15000, 45000, and 90000 (from left 
to right). 



Z = 10 within this paired configuration, whereas the 
corresponding unpaired strangelet has Z = 24. A sys- 
tematic study of the total charge of strangelets in this 
configuration will not be discussed here since this con- 
figuration is rather fragile with respect to the details of 
the single-particle spectra and thus difficult to realize for 
many different particle numbers. 

Let us therefore concentrate on the CFL-type solution, 
which exists for arbitrary particle numbers. We consider 
different mass numbers A from A — 108 to A ~ 90000, for 
one particular value of the coupling constant, HA^ — 2. 
In order to reduce the considerable numerical effort, we 
use for the large strangelets (starting from A = 15000) 
the condition (PO)) with /ie = (as a consequence, the 
quark numbers for each flavor and color are not inte- 
gers) instead of looking for the true energy minimum 
with respect to f3 decay. In addition, we do not mini- 
mize the energy with respect to the radius, but we sim- 
ply estimate the volume of the bag by dividing the mass 
number A by the baryon density pB bulk of infinite mat- 
ter. These two approximations are very accurate for such 
large strangelets. Already in the case oi A = 3000, the 
quark numbers and the radius are very well reproduced 
within these approximations: the full minimization re- 
sults in quark numbers N^r = Ndg ~ Ngt = 1052, 



ug 



N, 



dr 



Nub 



N,r = Nr,.h = Nsg = 974, and 



a radius R — 12.23 fm, while the approximations lead 
to Nur = 1051.8, Ndg = 1051.7, N,b = 1051.1, N^g = 
Ndr = 973.8, N^b = Nsr = 974.4, Ndb = N^g = 974.5, 
and R = 12.19 fm. Our results for the charge densities 
for A =108, 3000, 15000, 45000, and 90000 are shown in 
Fig.m 

Since all quarks are paired, we have equal numbers 
of up, down, and strange quarks such that the total 
charge of these strangelets is zero. The positive surface 
charge is mostly compensated by an excess of negative 
charge concentrated at around 1-3 fm below the surface. 
We stress that this concentration of negative charge in 
a thin layer is a consequence of pairing and the effect 
persists if Coulomb interaction is switched off. In fact. 



since also the strange quarks are paired, the "missing" 
strange-quark density at the surface must be compen- 
sated by an "overshooting" of the strange-quark density 
within a distance corresponding to the size of the Cooper 
pairs, i.e., the coherence length ^. Due to the strong gap, 
the coherence length is very small: Using the estimate 
^ ^ l/(7rA), one finds that it is of the same order as 
the Fermi wavelength and, strictly speaking, one might 
therefore question that mean-field results are quantita- 
tively correct [l^]. The smallness of explains why the 
compensation of the negative surface charge is mostly 
concentrated in such a thin layer at a small distance from 
the surface. 

Below this strongly negatively charged layer, the 
charge density stays negative but much smaller. Due to 
Coulomb interaction, which tries to push the charge to- 
wards the surface, this negative charge density decreases 
with increasing distance from the surface, especially for 
large strangelets. Actually, if Coulomb interaction is 
switched off, the remaining charge is distributed more 
or less homogeneously over the whole volume. 

The behaviour of the charge density far away from the 
surface in the presence of Coulomb interaction can eas- 
ily be interpreted in terms of Debye screening (similar 
considerations can be found in Ref. [l^ for the case of 
hadron-quark mixed phases): We know that in a uniform 
medium with Debye screening the Laplace equation for 
the Coulomb potential is replaced by 



V^-^]Ao = 0. 



(21) 



where A is the screening length, which can be obtained 
from the limit n°°(g° = 0,q ^ 0), where W{q) is the 
polarization tensor in the uniform system. This is equiv- 
alent to the expression [2^ 



dpch 
dpfc 



(22) 



Computing numerically this derivative within our model 
for the case of bulk CFL matter with S = 156.26 MeV 
and HA'^ = 2, we obtain A = 7.74 fm. 

Taking the Laplacian of Eq. (|2T|) . we see that the 
charge density obeys the analogous equation 



V'-^)Pc/.=0 



(23) 



In the case of half-infinite matter with a surface aX z = 
0, the solution of this equation shows that the charge 
density goes to zero as pch exp(z/A) if one goes away 
from the surface (z —oo). In the case of a sphere, the 
corresponding solution reads 



Pch oc 



sinh(r/A) 

T/a 



(24) 



Far away from the surface, the charge densities which 
we obtain are very well described by Eq. (|24p . To show 
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FIG. 9: Zoom into the part of Fig.[8]where the charge density 
behaves as given by Eq. (|24p . For a better visibihty, the charge 
densities have been divided by their respective values at r = 0. 
The thin dotted curve corresponds to Eq. (|24|l with A = 7.74 
fm. 



TABLE III; Fitted liquid-drop parameters for the CFL-type 
neutral strangelets [Z = 0). The surface tension a corre- 
sponding to the fitted value of as is also given. 



Bl/4 




as 


ac 




(MeV) 


(MeV) 


(MeV) 


(MeV) 


(MeV/fm^) 


156.26 


2 


107 


289 


11.9 


159.46 


2.25 


109 


297 


12.8 


163.46 


2.5 


112 


306 


13.9 



the rotated photon is massive (in fact, it is even heavier 



than the other gluons |27[). This means that in a large 
object, hke a strange star, all the negative charge will be 
concentrated within a layer of a thickness of at most a 
few tens of fm below the surface. However, before draw- 
ing any firm conclusion, one should study this problem 
in more detail. This will be left for future work. 



this, we display in Fig. [S] the same charge densities as in 
Fig. [51 but divided by their value at r = 0. Far away from 
the surface, all curves follow exactly Eq. (|24p with the 
value A = 7.74 fm calculated for bulk CFL matter. Near 
the surface, i.e., at distances which are of the order of a 
couple of Fermi wavelengths, there are strong deviations 
from this behavior due to Friedel-type oscillations [23 |. 
This is because Eq. (|2ip is not exact, but it is only valid 
in a uniform medium and in the long- wavelength limit. 

It is interesting to notice that the value of the De- 
bye screeing length we obtain is in reasonable agreement 
with the photon Debye mass calculated from perturba- 
tive QCD, which reads, for the CFL phase, m|, = 
l/A^ = 4 2i-8in2 g2jy^^2^(g^2-) ^5] {niD,-,^ dcuotes the 
Debye mass without giuon-photon mixing, see below). 
For typical values of the chemical potential this gives 
A - 10 fm. 

In principle, in color superconducting phases, the pho- 
ton can mix with one of the gluons. In the CFL phase, in 
bulk matter, one linear combination of photon and gluon 
stays massless. This means that at large distances d ^ £, 
the Debye screening for the "rotated" photon [2^, [23| 
does not work, since the Cooper pairs are neutral with re- 
spect to the rotated charge Q. Within the simple model 
we use for the moment, there are no gluons, such that 
the mixing cannot be studied. It could be taken into 
account, as mentioned at the end of Section [TTl by in- 
cluding the gluons in the same way as the photon, i.e., 
on the Hartree level. We expect that if we included the 
gluons in this way, we would find an even faster decrease 
of the charge if we go away from the surface, since in 
addition to the electromagnetic force we would have the 
color forces, which try to push the color charges to the 
surface, and in the CFL phase color neutrality goes hand 
in hand with electrical neutrality. Therefore, this is not 
in contradiction with the fact that the rotated photon is 
massless, but it is just a consequence of the fact that the 
combination of photon and gluon which is orthogonal to 



C. Liquid-drop type expansion 

The advantage of the present approach is that finite 
size effects are correctly implemented. For large num- 
bers of particles, this becomes, however, rather cumber- 
some and asymptotic expansions such as a liquid-drop 
type approach can be very useful. We will discuss here 
the determination of the parameters, such as the surface 
tension, of a liquid-drop type formula for the energy per 
baryon as a function of the baryon number A, including 
a surface and a curvature term, 

from our results. As in Section 1111 Bl we will restrict 
our discussion to the CFL-type solutions with Z = 0, 
such that we do not need to include a Coulomb term 

As explained after Eq. (f20|) . {E/A)i,uik should be the 
energy per baryon of infinite matter with — rather 
than that of (3 stable infinite matter. However, since 
we consider only the CFL-type solution, this distinc- 
tion is irrelevant. Hence, for our chosen parameter sets, 
we have {E/A)i,uik — 900 MeV. Since for the neutral 
strangelets the Coulomb interaction has only a negligible 
effect on the total energies (for example, in the case of 
the strangelets considered in Section fill Bi the Coulomb 
interaction changes the total energy per baryon by less 
than 5 keV) it will be neglected here in order to reduce 
the numerical effort. The result of the fitted coefficients 
as and ac for the different parameter sets are listed in 
Table lllll As an example, in order to show the accu- 
racy of the asymptotic expansion, we display in Fig. [TU] 
some results for the energy per baryon together with the 
liquid-drop formula, Eq. (|25p . The dashed line corre- 
sponds to the liquid-drop formula without the curvature 
term (ac = 0). From this figure it becomes clear that 
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FIG. 10: Energy per baryon as a function of baryon number 
for HA'^ = 2 and B^/" = 156.26 MeV. The exact results are 
indicated by the crosses, the fitted liquid-drop formula by the 
solid line. The dashed line corresponds to the liquid-drop 
formula without the curvature term. 



without color superconductivity [20| ■ On the other hand, 
the fact that the curvature term is very strong implies 
that the knowledge of the surface tension alone might not 
be sufficient in order to determine, e.g., the possibility of 
mixed phases, the size of droplets, etc. 

Before we conclude, let us comment on the physical 
meaning of the surface tension we obtain. In the MIT bag 
model, it is supposed that the energy needed to create a 
bag with volume V is simply given by BV. In principle 
one could imagine that there is an explicit dependence 
of the bag energy on, e.g., the surface or the curvature 
of the bag boundary. In Ref. this contribution to 
the surface tension was called "intrinsic surface tension" , 
(7/, and it was argued that it should be small. What 
we calculate here is the "dynamical surface tension", an, 
which has its origin in the change of the level density 
of the quarks inside the bag as a function of the bag 
geometry. 



the liquid-drop formula with curvature term works ex- 
treniely well, much better than in the case without pair- 
ing [28|. The reason is that shell effects are completely 
washed out because, contrary to the situation in ordinary 
nuclei, the pairing gap is much larger than the spacing 
between neighboring shells. Another interesting observa- 
tion is that the curvature term is very important, even 
for rather large mass numbers A. 

The coefficient as is closely related to a very interesting 
quan tity , namely the surface tension. As explained in 
Ref. 1231, the surface tension is obtained as 



where 



^5 



Es = E-A{^ 



(26) 



(27) 



bulk 



is the energy excess due to the surface and Rq is an ef- 
fective radius defined by 



PB bulk 



AttRI 



(28) 



which is actually very close to R for not too small 
strangclets. On the one hand, using the liquid-drop for- 
mula for E in Eq. one would obtain a surface 
tension which depends on A because of the curvature 
term. Therefore it is clear that one has to use Eq. (1^^ 
in the limit A — » oo, where the curvature term vanishes. 



(^SPslulk 

(36^)1/3 



(29) 



The corresponding numbers are given in the last column 
of Table IIIII They are of the same order of magnitude 
as the estimate a ~ (70 MeV)^ = 8.8 MeV/fm^ for SQM 



IV. SUMMARY 

In this paper we have investigated finite lumps of color 
superconducting SQM. To that end we have treated the 
MIT bag model, supplemented with a pairing interac- 
tion, in the framework of HFB theory. This allows us to 
correctly include finite size effects for pairing, too. The 
calculation is numerically rather involved, since in addi- 
tion to solving self-consistently the HFB equations, we 
have to determine the bag radius and the fractions of the 
different quark species by minimizing the total energy of 
the system. 

As expected from previous MIT bag-model studies, we 
find a suppression of the strange-quark densities at the 
surface, resulting in a positive surface charge. Our main 
result is that, in spite of this surface charge, the total 
charge of the CFL type solution is zero due to pairing, 
as in bulk matter. Most of the positive surface charge 
is compensated in a negatively charged layer situated at 
about 1-3 fm below the surface. The origin of this concen- 
tration of the negative charge is pairing: Since all quarks 
are paired, the positive surface charge must be compen- 
sated on a length scale corresponding to the coherence 
length. The remaining negative charge, which is neces- 
sary to compensate all of the positive surface charge, is 
situated below this layer. With increasing distance from 
the surface, the charge density decreases on a length scale 
of 8 fm, corresponding to the Debye screening length. 
This number will probably be strongly decreased if the 
gluons are included in a perturbative way similar to the 
photon. In any case, in the biggest part of a large object, 
such as a strange star, one finds vanishing charge density 
if one goes more than a few tens of fm away from the 
surface. It remains to be investigated in which way our 
results change the traditional picture of the surface of a 
strange star and the detectability of smaller strangelets 
in current experiments such as AMS-02 or LSSS |30 |. 

We have also compared our results for the energy per 
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baryon of finite strangelets with a liquid-drop like for- 
mula. We obtain a surface tension of the order of 12- 
14 MeV, in reasonable agreement with previous stud- 
ies where color superconductivity was not considered, 
and a strong curvature term which is crucial to repro- 
duce the correct energies up to baryon numbers of sev- 
eral thousands. An interesting result is that, in the 
presence of color superconductivity, the liquid-drop for- 
mula describes very accurately the total energies even for 
A < 100, at least for strangelets with even baryon num- 
ber. The reason is that, since the gap A is much larger 
than the spacing between the energy levels, shell effects 
are strongly suppressed. 
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Eq. Il]), translates into the following equation 

ffjnniR) ^ -gfj.niR) , (A5) 

or, explicitly. 



jliPfjnnR) = Sgn(Kn)./ "^^ jl'ipfjnnR) , (A6) 

where we number by n > the positive-energy (particle) 
states and by n < the negative-energy (antiparticle) 
states. In practice, we will keep only the states with 
positive eigenvalues and neglect the antiparticle contri- 
butions. The latter can approximately be absorbed into 
a redefintion of the coupling constant. 
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APPENDIX A: SPINORS IN A SPHERICAL 
CAVITY 

In this appendix we recall basic properties of free Dirac 
spinors in a spherical cavity (see, e.g., Ref. [Si). They 
can be written as 

where y are spinor spherical harmonics [31]. We have 
the following relations between the angular momentum 
quantum numbers 



= -0 + 2) 



^=3-\ ^'=J + l- (A2) 



For the solutions of the free Dirac equation, the functions 
/ and g are given as follows in terms of the spherical 

Bessel functions {CfjKn = \Jp}jKn + ™/) 

9fjK.n{r) —CfjK,njl{PfjK.nr) 



ffjnnir) =C/jKnSgn(Kn)./ ^-^^™ jl'ipfj^nr) , 



(A3) 



where the Cfj^n are normalisation coefficients which can 
be determined from the normalization 



drr^ / dn-i/j'' {r)il;{r) = 1 



(A4) 



The momenta p/jKn are obtained from the boundary con- 
dition. The boundary condition of the MIT bag model, 



APPENDIX B: HFB EQUATIONS 

In this appendix we will give some more details about 
the HFB equations. Their derivation is analogous to the 
derivation of the Dirac-Hartree-Bogoliubov equations in 
finite nuclei, which is given in Ref. hj. 

The HFB equations are derived from the Lagrangian 
by minimizing the energy in the mean field approxima- 
tion, i.e., linearizing the interaction under the assump- 
tion of nonzero expectation values for the condensates 
SAA'{x), Eq. ([7]). Due to the inhomogeneities of a finite 
system, the Green's functions become nondiagonal in mo- 
mentum. In the stationary case, it is convenient to work 
in r space for the spatial coordinates but to perform the 
Fourier transformation for the time variable. Then the 
Green's functions, 



with 



(Bl) 



(B2) 



take the following general form in Nambu-Gorkov space: 



G(r,r';c^) F(r,r';c^) 



Fir,r';u;) G(r,r';c^) 
J2 (^o.{r)\ 1 



q(£q>0) 



Va(r) I U! - Ea+irj 



(C/„(r'),T4(r')) 



E 

/3(e/3<0) 



Upir)\ 1 



{Uf,{r'),Vp{r')), (B3) 



where G, G and F, F are normal and anomalous Green's 
functions, respectively. The spinors Ua^p and Va^p corre- 
spond to the particle- and hole-like components, respec- 
tively. 
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The energy in mean-field approximation can now be 
written as 11911 



TT[{i-i ■ V - m)G{x,x+)] 



- / dSjTT[^ix,y)G{y,x+) - Aix,y)F{y,x+)] 



(B4) 

where the derivative in the first term acts only on x and 
not on and x^ means the four vector (a;*^ + t,x) in 
the limit t — > 0^. In our case, the normal and anomalous 
self-energies S and A are local and time- independent: 
■E{x, y) = eQAo(x)7°J(a:-y) and A(x, y) = A{x)S{x-y), 
and Eq. (|B4| can be reduced to Eq. (fT7)l . 

As mentioned in Section IIIB| the expectation values 
(like condensates, densities, etc.) which are needed for 
calculating self-consistently the self-energy E and the 
pairing field A can be expressed in terms of the U and 
V functions. To that end, it is sufficient to express them 
in terms of the Green's functions, e.g. 



SAA' 



-{i})T{x)TA\A'i^{x)) = iTrF(x,a;+)TAAA : 



(B5) 

which leads to Eq. Q. 

By minimizing the total energy with respect to the U 
and V functions, one obtains the HFB equations, [see 
Eq. ©]: 



(B6) 



with Wa = (UajVa)'^ and H being the matrix on the 
left-hand side of Eq. ([5]). 

For homogeneous infinite systems the matrix elements 
of Ti. are diagonal in momentum space and solutions to 
the HFB equations are known for many cases. For finite 
systems, in general, these equations are solved numeri- 
cally by diagonalizing the matrix Ti. in some conveniently 
chosen basis. Here, we are working in the basis which di- 
agonalizes the Dirac hamiltonian (i.e., hfc without the 
Coulomb potential), see Appendix 1X1 and the eigenvec- 
tors Ua{r) and Va{r) are developped within this basis. 

The matrix elements of the pairing fields A^(r) and 
of the Coulomb field Ao{r) are computed in the usual 
way. For illustration, we give here the explicit expression 
for the matrix elements of A2(r), which connects up and 
down quarks, in the basis described in Appendix [K\ 



(A2) 



jKnn' 



r<R 
R 



(r)A2(r)^ 



drr A2{r){gujnn{r)gdjK.n'{r) 



+ Iu]Kn{r)Id]Kn' [r)) . (B7) 

Note that, due to spherical symmetry, all matrices are 
diagonal in j and k and proportional to the unit matrix 
with respect to m. 



In spite of the spherical symmetry, the matrix to be 
diagonalized is still huge, limiting the baryon number 
which can be calculated with reasonable computational 
effort. It is therefore important to reduce the size of 
the actual matrix to be diagonalized. By means of an 
orthogonal transformation 



n = sns^, w = sw, ss^ = i 



(B8) 



in color, flavor, and Nambu-Gorkov space, the mateix 
can actually be block-diagonahzed (see, e.g., Ref. js j |) 
containing seven blocks. Six of them, Ti.B....Gi ^-re 2 x 
2 matrices in Nambu-Gorkov space, describing mutual 
pairing of two particles, such as, e.g., red down quarks 
(dr) with green up quarks (ug): 



Hi 




(B9) 



where hfc is the single particle Hamiltonian for flavor / 
and color c. The second and third 2x2 blocks are 




Ti-D — 




(BIO) 



Since we have in addition the pairwise relations 'He,f,g = 
—Hb,c,d, only three of the six 2x2 blocks have to be 
diagonalized in practice. The seventh block. Ha, is 6 x 6 
in Nambu-Gorkov space and describes pairing between 
red up, green down and blue strange quarks 



Ha 
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-hdg 





VAs 
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(Bll) 



APPENDIX C: CUTOFF FOR THE GAP 
EQUATION 



As mentioned in Section [II Bl the divergent gap equa- 
tion is regularized with the help of a smooth cutoff func- 
tion 



/(P/A) = 



1 



1 -I- ci exp(c2a(p/A - 1)) 



(CI) 



where ci = \/2 - 1, C2 = 1/(4 - 2^2), and a = 22.58 
have been chosen such that (p/A) approximates the 
cutoff function g{p/A) used in Ref. [33|, but our function 
has the advantage to fall off more rapidly at very high 
momenta, which allows us to truncate the basis at a lower 
energy. 

This function is used as a form factor multiplying each 
of the four legs of the four-point vertex. In practice, 
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this means that the form factor is used in two places: 
First, when calculating s^yi(r), and second, when calcu- 
lating the matrix elements of A^(r) in the basis of the 
spinors defined in Appendix \X\ It should be noted that 
the diagonalization of the HFB matrix does not directly 
provide us with the eigenfunctions Uair) and VQ(r), but 
with their respective expansion coefficients in the basis 



of the spinors defined in Appendix [X] When calculating 
SAAi^) according to Eq. the coefficients have to be 
multiplied with the corresponding basis functions, and in 
this step the factor /(p/j^n/A) is attached to each basis 
function. Second, when calculating the matrix elements 
of the gap Aa, we again attach a factor fipfj^n/-^) to 
each basis function. 
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